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The spreading of infectious diseases with and without immunization of individuals can be modeled 
by stochastic processes that exhibit a transition between an active phase of epidemic spreading and 
an absorbing phase, where the disease dies out. In nature, however, the transmitted pathogen may 
also mutate, weakening the effect of immunization. In order to study the influence of mutations, 
we introduce a model that mimics epidemic spreading with immunization and mutations. The 
model exhibits a line of continuous phase transitions and includes the general epidemic process 
(GEP) and directed percolation (DP) as special cases. Restricting to perfect immunization in two 
spatial dimensions we analyze the phase diagram and study the scaling behavior along the phase 
transition line as well as in the vicinity of the GEP point. We show that mutations lead generically 
to a crossover from the GEP to DP. Using standard scaling arguments we also predict the form 
of the phase transition line close to the GEP point. It turns out that the protection gained by 
immunization is vitally decreased by the occurrence of mutations. 
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I. INTRODUCTION 

The modeling of epidemic spreading is a fascinating 
subject both in theoretical biology and statistical physics 
far from equilibrium 0, @ • A possible approach is the 
study of stochastic models that mimic the competition 
of infectious spreading and recovery by certain proba- 
bilistic rules. Depending on the rates for infection and 
recovery the disease may either spread over the popula- 
tion or disappear after some time. 

The simplest models for epidemic spreading assume 
that the individuals live on the sites of a (^-dimensional 
lattice. At a given time each individual can be either 
infected or healthy. The system evolves according to cer- 
tain probabilistic rules that resemble infection of nearest 
neighbors and spontaneous recovery. If the susceptibility 
to infections is sufficiently large the epidemic will spread 
while for low susceptibilities it will die out. 

In most models it is assumed that the contagious dis- 
ease is transmitted exclusively by direct contact. This 
means that the disease, once extinct, cannot appear 
again. The system is then trapped in a fully recovered 
state, which can be reached but not be left. Such a state, 
where the system is dynamically trapped, is called ab- 
sorbing. The two regimes of spreading and extinction are 
usually separated by a so-called absorbing phase transi- 
tion. In the supercritical regime, where infections dom- 
inate, the epidemic may spread over the entire system, 
reaching a fluctuating stationary state with a certain av- 
erage density of infected individuals. In the subcritical 
regime, where recovery dominates, the system eventually 
reaches the fully recovered absorbing state. 
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Close to the transition the temporal evolution of the 
spreading process is characterized by large-scale fluctua- 
tions. Theoretical interest in epidemic spreading stems 
from the fact that this type of critical behavior is univer- 
sal, i.e., it does not depend on the details of the model un- 
der consideration. The classification of all possible tran- 
sitions from fluctuating phases into absorbing states is 
currently one of the major goals of nonequilibrium sta- 
tistical physics 

In epidemic models without immunization, where the 
pathogen (e.g. a virus) is transmitted to nearest neigh- 
bors, the critical behavior at the transition generically 
belongs to the universality class of directed percolation 
(DP) jfj. The DP class is very robust and plays an im- 
portant role not only in epidemiology but also in various 
other fields such as Reggeon field theory Q , interface de- 
pinning , population growth § , catalytic reactions [tj , 
and flowing sand llOl (for a review of possible experimen- 
tal applications see |11|). 

As a next step towards a more realistic description of 
epidemic spreading one can take the effect of immuniza- 
tion or weakening by infections into account 0, IT^L IT3 , 
HI3. The simplest way of implementing immuniza- 
tion is to change the initial susceptibility of an individual 
after the first infection and keep it constant thereafter. 
Generally such a process is controlled by two parameters, 
namely, a probability for first infections and a reinfection 
probability. In the case of perfect immunization, where 
each individual can be infected only once, one obtains 
the so-called 'general epidemic process' (GEP) 0,0, li- 
lt differs from DP in so far as the disease can only spread 
in those parts of the system which have not been infected 
before. Thus, starting with a single infected site in a non- 
immune environment, the disease typically propagates in 
form of a solitary wave, leaving a region of immune sites 
behind. Depending on the infection rate this wave may 
either spread over the entire system or stop before. The 
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transition between infinite and finite spreading is a criti- 
cal phenomenon which in this case belongs to the univer- 
sality class of dynamical isotropic percolation . Note 
that unlike DP models, a GEP running on a finite sys- 
tem has no fluctuating active state, instead the process 
terminates when it reaches the boundaries. 

In nature, however, immunization is a much more com- 
plex phenomenon. For example, the protection by immu- 
nization may abate as time proceeds. Even more impor- 
tantly, the strategy of immunization competes with the 
ability of the contagious pathogen to mutate so that it 
can no longer be recognized by the immune system of 
previously infected individuals, weakening the effect of 
immunization. The aim of the present work is to intro- 
duce and study a simple model which mimics epidemic 
spreading with immunization and mutation. To this end 
we generalize the model described in [H by includ- 
ing mutations as well as a mechanism for the competi- 
tion between different species of pathogens. The model 
is controlled by three parameters, namely, a first infec- 
tion probability, a reinfection probability controlling the 
effect of immunization, and a probability for spontaneous 
mutations. Moreover, the model is defined in such a way 
that it includes DP and GEP as special cases. 

The article is structured as follows. In Sec.[n]we first 
define the model. Restricting our analysis to the case 
of perfect immunization in two spatial dimensions, we 
discuss the phase diagram and the qualitative behavior 
of the model in Sec. IIIII In Sec. IIVI the critical behav- 
ior at the phase transition line is studied in detail while 
crossover phenomena in the vicinity of the GEP point are 
investigated in Sec. |VJ The article ends with concluding 
remarks in Sec. IVII 



II. SIMULATION MODEL 

Our model is meant to describe the spreading of an in- 
fectious disease that evolves as follows. Individuals can 
be healthy or infected with a certain pathogen. During 
their illness infected individuals may infect or reinfect 
neighboring individuals with certain probabilities. More- 
over, there is a probability that a pathogen mutates dur- 
ing transmission. Because of the enormous number of 
possible mutations one usually obtains an entirely new 
type of pathogen which has not been involved before. To 
simplify the model we also assume that each individual 
can be infected at a given time by no more than a single 
type of pathogen. If the individual is exposed simulta- 
neously to several competing pathogens, one of them is 
randomly selected. 

In more technical terms the model is defined as follows. 
Individuals live on the sites of a <i-dimensional simple 
cubic lattice. Each site may be either healthy (inactive) 
or infected (active) by a pathogen of type n, where n is 
a positive integer. Moreover, each site keeps track of all 
species of pathogens by which it has been infected in the 



past. Therefore, the state of a site is characterized by a 
number n together with a dynamically generated list of 
all previous types of infections. 

The model evolves in time by synchronous updates, 
i.e., in each time step the whole lattice is updated in 
parallel as follows: 

1. Spreading: Each infected individual at time t 
transmits its pathogen n to its 2d nearest neigh- 
bors (target sites) at time t + 1. 

2. Susceptibility: The transmitted pathogen 
reaches the target site at time t + 1 with proba- 
bility p if this site has never been infected before 
by this species, otherwise it reaches the target site 
with probability q. 

3. Competition: If a target site is exposed to several 
transmitted pathogens, one of them is randomly 
selected with equal weight. 

4. Mutation and Infection: Before infecting the 
target site the selected pathogen mutates with 
probability A, replacing n by a new integer num- 
ber (drawn from a global counter) which has not 
been used before. 

5. Immunization: In case of a first infection the type 
of pathogen is added to the list of species against 
which the site will be immune in the future. 

6. Recovery: All sites that have been active at time 
t recover until time t + 1 unless they are not again 
infected during the update from t to t + 1. Hence 
the time of illness is a single time step. 



III. PHASE DIAGRAM 

In what follows we restrict our analysis to the special 
case of perfect immunization q = 0. In this case the 
model is controlled by only two parameters, namely, the 
probability of first infections p and the probability of mu- 
tations A. Moreover, we restrict ourselves to the case of 
d = 2 spatial dimensions. The corresponding phase di- 
agram is shown in Fig. ^ It comprises an active phase, 
where the epidemic spreads, and an inactive phase, where 
the disease dies out so that the system eventually enters 
the fully recovered absorbing state. Both phases are sep- 
arated by a curved phase transition line. 

We first note that the endpoints of the phase tran- 
sition line correspond to well-known special cases. On 
the one hand, for A = the model reduces to the GEP 
on a square lattice @, provided that only one type of 
pathogen is involved. In this case the critical value of p 
is exactly given by p c = 1/2 16]. On the other hand, for 
A = 1 all transmitted pathogens mutate, i.e., the target 
sites are always infected with a new species so that im- 
munization has no influence. It is easy to see that in this 



3 




X 



FIG. 1: Phase diagram of the model. The active and the 
inactive phases are separated by a curved line of continuous 
phase transitions, connecting the points of critical GEP and 
critical DP (see text). Circles mark the numerically deter- 
mined critical points. The inset shows the phase transition 
line for small A. 

case the model reduces to directed bond percolation on 
a square lattice with p c ~ 0.28734 ^t|. The other points 
on the phase transition line in Fig. ^ were determined 
numerically using seed-simulations (see Sec. IIV A"|) . 

As already discussed in the Introduction, mutations 
generally weaken the effect of immunization. This ex- 
plains why the critical value of p decreases monotonically 
with increasing mutation probability A. Fig. ^ clearly 
shows that introducing mutations in a GEP has a strong 
influence on p c , especially if A is very small. Since the 
phase transition line approaches the GEP point with infi- 
nite slope (see inset in Fig.^l a tiny increase of A reduces 
the corresponding critical value of p dramatically. This 
indicates that mutations are a relevant perturbation and 
thus the spreading behavior of the model for A > is 
expected to differ from that of the GEP. On the other 
hand, for larger values A > 0.1 the critical value of p 
decreases only moderately with increasing A, indicating 
that the behavior of the model in this region is essentially 
the same as for A = 1, where the transition is known to 
belong to DP. 

Fig. |21 shows snapshots of simulations at the critical 
point for different times and various values of A. Infected 
individuals are represented by black dots. Moreover, in- 
dividuals which are immune against at least one active 
type of pathogen are marked by gray dots. If a site is 
immune solely against pathogens that already became 
extinct the gray dot is removed. The snapshots shown in 
this figure suggest the following qualitative behavior: 

a) For A = the spreading process creates a growing 
cluster of immune sites with infected individuals 
located at the edges. As usual in the GEP, the re- 
gion of immune sites is compact in the active phase 




b) 



c) 



d) 



FIG. 2: Snapshots of seed-simulations at the phase transition. 
The black square marks the position of the seed. Each row 
shows a simulation for a fixed value of A and times t = 30, 
t — 100 and t = 225 (from left to right). From a) to d) values 
of A are (GEP), 0.003125, 0.05 and 1 (DP). Black (gray) 
dots denote active (immune) individuals (details in the text). 



while it is fractal at the phase transition. 

b) For a small value of A the disease first behaves like a 
GEP (for t < 100) while for larger times the spread- 
ing behavior changes and clearly differs from the 
GEP. Although there still is a region of immune 
individuals it does no longer provide efficient pro- 
tection against infections since it is reinvaded by 
mutated pathogens. 

c) Increasing A further the process changes its appear- 
ance already at an early stage. There are only small 
patches of immune sites and the process reminds 
more of DP than GEP. 

d) Finally for A = 1 every transmitted pathogen mu- 
tates into a new one. In this case immunization has 
no influence and the process reduces to DP. Since 
the disease is not driven away from immunized re- 
gions the spatial extension of an epidemic at the 
phase transition grows much slower compared to a 
critical GEP. 
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Based on these phenomenological observations we expect 
the model to behave initially in the same way as a GEP. 
After a certain time mutations become relevant, allowing 
former immune areas to be reinvaded. Especially close to 
the transition the process survives long enough to reach 
this crossover time. The visual appearance of the process 
is then increasingly similar to that of a DP process. The 
time it takes to observe the crossover from GEP to DP 
grows with decreasing A and eventually diverges in the 
limit A -> 0. 



< A < 1 we note that this would be a non-trivial result 
since the so-called DP-conjecture does not apply in the 
present case. The DP-conjecture 0] states that phase 
transitions in two-state systems with a reachable absorb- 
ing state and short-range interactions belong to DP, pro- 
vided that memory effects, non-conventional symmetries, 
and quenched disorder are absent. Contrarily the present 
model has many absorbing states and memorizes previ- 
ous infections over a long time. 



IV. CRITICAL BEHAVIOR ALONG THE 
PHASE TRANSITION LINE 



Phase transitions into absorbing states are usually 
characterized by simple scaling laws. In models for epi- 
demic spreading an important quantity is the probabil- 
ity P s (t) that an epidemic starting from a single infected 
seed in a healthy and non-immune environment survives 
at least until time t. 

Let us first recall the standard scaling laws for the sur- 
vival probability. Let A p = p — p c denote the distance 
from criticality. If |A p | -c 1 the survival probability 

t~ b until a certain 
Ap| _y n is reached from where on it ei- 

8v\ 



first decays as a power law P s {t) 
timescale £| 

ther saturates at P s (oo) ~ Ap ^ 11 for A p > or decays 
exponentially for A p < 0, reaching a spatial extension 
£j_ ~ lApl - ^. Assuming scaling invariance this behav- 
ior can be described in terms of a scaling form 



(i) 



where $(0 is a scaling function with the asymptotic be- 
havior 
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such that the time dependence drops out for £ — > +oo. 
The scaling function and the triplet of critical expo- 
nents (i/\\,i/±,5) are believed to characterize the univer- 
sality class of the phase transition under consideration. 
We note that this scaling form is known to be valid both 
for GEP and DP, although with different sets of critical 
exponents 



(1.506, 4/3, 0.092) for GEP 
(1.295, 0.734, 0.451) for DP. 



(3) 



For the density p in DP a similar scaling form as in Eq. Q 
is valid (with the same exponents as in Eq. J3J for d = 2). 

We now analyze the critical behavior of the model, as- 
suming that the scaling form Q is valid everywhere in 
the vicinity of the phase transition line. Although the 
qualitative discussion in Sec. |U suggests DP behavior for 



A. Seed simulations 

Seed simulations start with a single infected site at 
the origin in a non- immune environment. Each run is 
stopped when it dies out or reaches a preset maximum 
time. We average over many runs with different real- 
izations of randomness. In order to eliminate finite size 
effects the lattice is always chosen large enough so that 
the epidemic never reaches its boundaries. As usual in 
this type of simulation we measure the survival probabil- 
ity P s (t), the number of active sites N(t) averaged over 
all runs, and the mean square spreading from the origin 
R 2 (t) averaged over all active sites in surviving runs |2Cj . 
The scaling form Q implies that these quantities vary 
at criticality algebraically as 



P B (t) ~ t 



-5 



N(t) ~ t 6 , R 2 (t) ~ t 



2/z 



(4) 



where z = v\\/v±_ and 9 = d/z — 26 — 1 for GEP and 
9 = d/z - 25 for DP [2l[. 

Fig. |31 shows our results of seed-simulations along the 
phase transition line. We note that these simulations are 
numerically challenging since for each site a list of infec- 
tions in the past has to be created and continuously up- 
dated. Therefore, numerical simulations are mainly lim- 
ited by the available memory and restricted the present 
study to a maximal temporal rang e of about 10 4 Monte 
Carlo updates. 

In order to determine the critical threshold p c {\) we 
kept A fixed and varied p until P s (t) displayed the ex- 
pected slope of DP in a log-log plot (dashed lines in 
Fig. [3J . Although this procedure is in favor of DP scaling 
the mere fact that it works consistently for all quantities 
in Eq. confirms that the transition does belong to 
DP for any < A < 1 while GEP scaling can be ruled 
out. The data also shows the expected crossover. For 
small A the curves roughly display the slope expected for 
GEP (dotted lines) before they cross over to DP, con- 
firming the crossover scenario discussed in the previous 
sections. Thus the introduction of mutations in a GEP is 
a relevant perturbation in the sense that it changes the 
asymptotic critical behavior of the model. At the phase 
transition the interplay between immunization and mu- 
tations drives the system towards DP. 

We note that the involvement of different species of 
pathogens in a spreading process with mutations allows 
one to introduce the number of active species N sp (t) as 
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an additional order parameter. As shown in panel d) of 
Fig- El the number of active species at criticality increases 
in the same way as the number of infected individuals. 
Thus their quotient tends to a A-dependent constant 

c(A) = lim jV/JV sp . (5) 

t— >oo 



B. Full-lattice simulations 

In this type of simulation we use a finite system of size 
L x L with periodic boundary conditions. The initial 
configuration is a fully occupied lattice where all indi- 
viduals are infected by different types of pathogens. (It 
would be also possible to occupy all sites with the same 
type of pathogen, but then the process would immedi- 
ately be trapped in the absorbing state.) We measure 
the density of active sites p{t) and the density of active 
species p sp (t). In the case of directed bond percolation 
A = 1 it is known that the density of active sites decays 
as p(t) ~ t~ s . Exemplarily we performed a full- lattice 



simulation for A = 0.5 at the critical point, confirming 
this type of decay with S = 0.451. Moreover, like in seed 
simulations, the density of active species also decays as 

The three exponents S, 6, z in Eq. that govern 
the spreading behavior at the transition depend only on 
two of the three independent critical exponents that are 
needed to characterize the universality class of DP (see 
Eq. ©)■ In order to check the value of the third indepen- 
dent exponent we performed off-critical full-lattice simu- 
lations for A = 0.5 and different values of < A p -C 1. 
As is shown in Fig. fusing the critical exponents 6 and i>» 
of DP one obtains a reasonable data collapse of the den- 
sity p(t). This supports the claim that along the phase 
transition line for A > indeed all three exponents gov- 
erning the scaling behavior of the spreading process are 
that of DP. 
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FIG. 4: Data collapse of p(t) for A = 0.5 using the critical 
exponents of DP, i.e., 8 — 0.451 and u« = 1.295. 



CRITICAL BEHAVIOR IN THE VICINITY 
OF THE GEP POINT 



In this section we investigate the influence of muta- 
tions in the vicinity of the GEP point in order to address 
two questions, namely, how does the system cross over 
from GEP to DP and why does the phase transition line 
terminate in the GEP point with an infinite slope. 



A. Decay at the GEP point 

Already at the GEP point our model exhibits a new 
feature, namely, the competition of different types of in- 
fections in full-lattice simulations (see Sec. lIVB|) . Start- 
ing the simulation in a finite system, where each site is 
initially infected by a different kind of pathogen, one ob- 
serves a coarsening process of competing species, leading 
to a slow decay of the density of active sites and active 
species. As shown in Fig [SI the decay of these quanti- 
ties suggests possibly asymptotic power laws, although 
the data display a considerable curvature in both cases. 
Assuming asymptotic power laws 



P (t)~t- a A»(t)~r a (6) 

at the GEP point we extrapolate the effective exponents 
a(t),a(t) visually for t — > oo (see insets of Fig. 0, ob- 
taining the estimates 



a = 0.66(2) , 



1.33(4) 



(7) 



suggesting that a = 2a. 

Regarding the limited accuracy of our numerical sim- 
ulations the conjecture of an asymptotic algebraic de- 
cay has to be taken with care. However, indirect sup- 
port comes from the one-dimensional case. Here the 
GEP transition is shifted to p c = 1 and the dynamics 
of competing species reduces to a ballistic coalescence 
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FIG. 5: Simulations at the GEP point p = 1/2, A = starting 
with a fully occupied lattice, where each site carries a different 
kind of infection. Although the lines are slightly curved, the 
decay of p and p sp may suggest a possible power-law behavior 
in the limit t — > oo. As shown in the insets, the correspond- 
ing exponents can be extrapolated visually by plotting the 
effective slope against 1/y/t (see text). 



process [23, for which asymptotic power laws could be 
derived exactly. 



B. A-controlled transition at the GEP point 

Let us now turn to the critical behavior in the vicin- 
ity of the GEP point. The GEP point in Fig. □ can be 
approached either vertically by varying p or horizontally 
by varying A. The critical behavior in vertical direction 
has been studied in detail in fiefs. H, 0, ^| and can be 
described in terms of the scaling form Contrarily, 
moving in horizontal direction by varying A and keeping 
p = 1/2 fixed one encounters mutations as a new feature, 
leading to a non-trivial fluctuating active state. 

As usual in critical phenomena, the additional con- 
trol parameter A is associated with a novel critical expo- 
nent Like vu this exponent is defined in such a way 



FIG. 6: Data collapse of the survival probability. 



that 



f||~A- 



(8) 



is the correlation time in the stationary state for p = 1/2 
and < A < 1. Similarly, the corresponding spatial 
correlation length is expected to scale as 



A" 



(9) 



where = p\\/z. To verify this conjecture we measured 
the survival probability P s (t) in seed simulations at p = 
1/2 for two different values of A. According to standard 
scaling theory this quantity should obey the scaling form 



P s (t) st _J *(At 1 /' i «), 



(10) 



where S ~ 0.092 is the density decay exponent of the 
GEP. Eq. I|10|) implies that the survival probability even- 
tually saturates at a value P s (oo) ~ A^n . Using this 
scaling form we collapse the data sets in Fig. El obtain- 
ing the estimate 



= 0.63(3) . 



(11) 



Turning to full-lattice simulations we can combine this 
scaling law with the results of the previous subsection, 
leading to the conjecture that p(t) and p sp (t) obey the 
scaling forms 



p(t) = t-^At 1 /"") 

Psp (t) = t- & n(M^^), 



(12) 
(13) 



where 0(C) and 0(C) are scaling functions with the 
asymptotic behavior 



0(0 = 




(14) 



Again this implies that eventually the densities reach sta- 
tionary values p(oo) ~ A QAI n and p sp (oo) ~ A QM n. Al- 
though there are deviations for small t these scaling forms 
lead to reasonable data collapses, as shown in Fig. [7| 
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FIG. 7: Data collapse of p and p sp based on the scaling forms 
112II . 113H using the exponents of Eqs. P. 1111 . 



We now suggest an explanation for the numerically de- 
termined value of /iii = 0.63. Initially the process behaves 
as a critical GEP until mutations become relevant at a 
typical time C||- Our argument is based on the assump- 
tion that £|| scales in the same way as the typical time 
at which the first mutation occurs. With the mutation 
probability A one needs on average A -1 infections until 
the first mutation occurs. As the process initially be- 
haves as a critical GEP the number of infections grows 
as / dtN(t) ~ t e+1 . Hence we are led to 



C|| - A" 



(15) 



with 9 = d/z — 2(5 — 1 = 0.587. This implies the scaling 
relation 



1 1.587 



= 0.630 



(16) 



which is in perfect agreement with the numerical estima- 
tion in Eq. {TTJ . 



C. Curvature of transition line at the GEP point 

So far the numerical analysis suggests that in the vicin- 
ity of the GEP point the epidemic process with mutations 
is invariant under scaling transformations of the form 



x — > x' = bx 

t t' = b z t 

A p -► A p ' = ft" 1 /"-'- A p 

A -> A' = 6 _1/,/1 - L A , 



(17) 
(18) 
(19) 
(20) 



where 6 is a scaling factor, A p = p— 1/2, (5 and 2 are the 
critical exponents of GEP, and fj,± — fi\\/z. In addition, 
the order parameters have to be rescaled appropriately. 
In seed simulations this means that 

P s (t)^P' B (t') = b- 5z P s (bH), (21) 

leading to the combined scaling form 

P s (t, A p , A) = r 5 $(A p t 1 < /l 'ii , At 1 /"" ) (22) 

in the vicinity of the GEP point. Similarly, in full-lattice 
simulations one would have to rescale 

p(t) - p\t') = b- az p(b z t) (23) 
Ps P (t) - P ' ap (f) = b- &z Psp (b z t) (24) 

provided that the conjecture of asymptotic power-law be- 
havior in Eq. © is correct. This would then lead to the 
scaling forms 

p(t,A p ,X) = t~ Q 0(Ap£ 1 / l/ ii , At lAl " ) (25) 
p S p(i,A p ,A) = t-^Apt^^Xt 1 /^). (26) 

As usual in the theory of critical phenomena, the phase 
transition line itself has to be invariant under scaling 
transformations. Comparing Eqs. (|19fl and (|2(Jfl we are 
led to the conclusion that the form of the transition line 
for small values of A is given by 



where 



A r 
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V\ V\\ 



A 7 . 



0.63 
1.506 



= 0.42. 



(27) 



(28) 



Since 



9A D 



A 



-0.58 



-g^- ■ - /\ the phase transition line indeed ter- 
minates at the GEP point with an infinite slope. 

In order to confirm the relations (|27l) . 128|) we plotted 
|A p | versus A -C 1 in a double-logarithmic representation 
in Fig.|S| The local slope of this curve leads to the effectve 
exponent (inset in Fig. |SJ which can be extrapolated and 
leads to the estimation 7 w 0.41(3), in agreement with 
the prediction in Eq. I|28|) . 
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FIG. 8: Double-logarithmic plot of the phase transition line 
|A P (A)| for AC 1. The local slope of this curve leads to the 
effective exponent which is shown in the inset together with 
its error bars. The extrapolation of the effective exponent for 
a/X — > leads to an estimation of the value of 7 in Eq. I27H 
which is in agreement with the prediction of Eq. 12811 . 



VI. DISCUSSION AND CONCLUSION 

In this paper we have introduced a minimal model for 
epidemic spreading with immunization and mutations. 
Apart from the probabilities for first infections and rein- 
fections p and q the model is controlled by a probability 
A that a transmitted pathogen mutates, creating a new 
pathogen which was not involved before. The model in- 
cludes the GEP (A = 0, q = 0) and DP (A = 1 or q = p) 
as special cases. 

Restricting the analysis to the case of perfect immu- 
nization and two spatial dimensions we have shown that 
the transition between survival and extinction in our 
model belongs to the universality class of DP everywhere 
along the phase transition line except for the point where 
the model reduces to critical GEP. In the vicinity of 
the GEP point even a small mutation probability drives 
the system away from criticality into a fluctuating active 
state. The crossover from GEP to DP can be described in 
terms of a suitable scaling theory, which involves a new 
exponent jiu. This crossover exponent also determines 
the form of the phase transition line in the vicinity of the 
GEP point. We suggested an explanation for the value 
of which turns out to be in perfect agreement with 
the numerical analysis. 

Although the model presented here is highly idealized 
(individuals on a square lattice, homogeneous infection 
probabilities, nearest-neighbor infections etc.) there is 
an important conclusion to be drawn regarding realistic 
spreading of epidemics in nature. As in the model, real- 
istic epidemic spreading starts at a certain threshold de- 
termined by various parameters such as the average sus- 
ceptibility, the interaction frequency, and the degree of 
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immunization and/or vaccination. Mutations weaken the 
effect of immunization, thereby decreasing this thresh- 
old. An important message of our paper is that for a 
population which is mainly stabilized by immunization 
and/or vaccination this threshold varies nonlinearly with 
the mutation rate, in the present case roughly as the 
square root of A. Thus even a small rate of mutations 
can significantly weaken the stability of a population at 
the onset of epidemic spreading. 

As a possible extension of the present study it would 
be interesting to investigate whether these properties can 
also be observed in epidemic processes with long-range 



infections 23] if mutations are introduced. Moreover, it 
would be interesting to stud y th e surface critical behavior 
at the system's boundaries |24| . 
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